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Abstract. We consider the counting statistics of electron transport through a double 
quantum dot with special emphasis on the dephasing induced by a nearby charge 
detector. The double dot is embedded in a dissipative enviroment, and the presence of 
electrons on the double dot is detected with a nearby quantum point contact. Charge 
transport through the double dot is governed by a non-Markovian generalized master 
equation. We describe how the cumulants of the current can be obtained for such 
problems, and investigate the difference between the dephasing mechanisms induced 
by the quantum point contact and the coupling to the external heat bath. Finally, we 
consider various open questions of relevance to future research. 
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1. Introduction 

The study of random fluctuations has a relevant role in many branches of physics [H 
[21 [3l H]. Close to equilibrium, fluctuations are intimately connected with dissipative 
relaxation mechanisms according to the fluctuation-dissipation theorem, independently 
of the physical origin of the fluctuations, classical or quantum-mechanical [3 E]. In 
contrast, far from equilibrium fluctuating quantities provide a unique insight into the 
internal properties of the system under consideration [7]. An immediate example is the 
evaluation of the quasi-particle charge of carriers through measurements of the current 
cumulants [8]. 

Many important phenomena can be characterized in terms of counted, elementary 
entities. The concept of particles, in the quantum realm, naturally defines what 
quantities should be counted. From this point of view one recognizes that counting 
problems constitute a quite general framework in which many different dynamical 
processes can be interpreted, also in the presence of complex quantum physics. The 
first application of the counting approach in quantum physics came from photon 
counting experiments, where the concept of full counting statistics (FCS) was originally 
developed [9]. Recently, this concept has attracted intensive theoretical [1] and 
experimental [lOl |Tll [12] attention within the field of electron transport. In the context 
of mesoscopic transport, FCS was introduced in order to characterize the noise properties 
of nanodevices [13]. Later, it was demonstrated also to be a sensitive diagnostic tool for 
detecting quantum-mechanical coherence, entanglement, disorder, and dissipation [1]. 

Mathematically, FCS encodes the complete knowledge of the probability 
distribution P{n, t) of the number n of transmitted entities during the measurement 
time t or, equivalently, of all corresponding cumulants. The study of counting statistics 
for stochastic processes is generally of broad relevance for a wide class of problems. 
For example, non-zero higher order cumulants describe non-Gaussian behavior and 
contain information about rare events, whose study has become an important topic 
within non-equilibrium statistics in physics, chemistry, and biology |6l [HI [15] . Within 
the framework of master equations some important results were recently obtained. 
Bagrets and Nazarov [16] have shown that the cumulant generating function (CGF) 
corresponding to a Markovian master equation is determined by the dominating 
eigenvalue of the rate matrix, when counting fields are appropriately included. Some of 
us have shown that it in principle is possible to calculate arbitrary orders of cumulants 
using perturbation theory in the counting field, rather than solving the full eigenvalue 
problem for the dominating eigenvalue [T7]. For non-Markovian systems described by 
generalized master equations (GME), we have shown that the CGF scales linearly with 
time [T8| [19], as in the case of Markovian processes, if the memory kernel has no 
power-law tails. Moreover, the CGF can be calculated using a so-called non-Markovian 
expansion [18] . 

Recently, we developed a method which unifies and extends these earlier approaches 
to FCS within a GME formulation [19]. Due to their intrinsic analytic structures, the 
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previous approaches were in practice limited to systems with only a few states [TBI [TS] . 
or only the first few current cumulants could be addressed [17]. In contrast, this recent 
advancement enables studies of a much larger class of problems, including the evaluation 
of zero-frequency current cumulants of very high orders for non-Markovian systems with 
many states p2]. We also showed how the method allows calculations of the finite- 
frequency current noise for non-Markovian transport processes [191 120]. A detailed 
account of these techniques will be given elsewhere [21]. In this paper we mainly want 
to restate the essential findings and address the open questions of relevance in future 
research. 

In order to demonstrate the applicability of our methods, we consider the current 
fluctuations of charge transport through a coherently coupled double quantum dot, 
a charge qubit. We consider the effects of a nearby quantum point contact (QPC) 
charge detector and the coupling to a dissipative phonon bath. If the QPC barrier 
height is modulated by electrons on the double quantum dot, the fluctuations of the 
current through the QPC monitors the dynamics of the qubit. This introduces a 
qubit dephasing mechanism. As we shall see, current fluctuations can be useful for 
extracting information about the internal dynamics of the double dot system. We 
concentrate on the transition between coherent and incoherent transport through the 
qubit, demonstrating the sensitivity of the cumulants to this transition. 

The structure of the paper is as follows: In Section [2] we summarize the general 
concepts of our method, while clearly identifying the essential steps to obtain the PCS 
or the current cumulants for a system governed by a non-Markovian GME. We briefly 
sketch the derivation of the expressions for the first few current cumulants (current, 
noise, and skewness) used in the following section. In Section [3] we describe the model 
of a dissipative qubit with a nearby QPC charge detector and show how the current 
cumulants yield information about the dynamics of the qubit. In Section H] we discuss 
various open questions and give an outlook for future research. 

2. Non-Markovian GME 

For many nanoscopic systems it is convenient to consider the evolution of only a few 
degrees of freedom. The system evolution is then captured by the dynamic equation for 
the reduced density matrix corresponding to these degrees of freedom. This equation 
should contain the effects of all external forces driving the system and, at the same time, 
the effective dynamics due to the degrees of freedom that have been traced out. The 
dynamics of the reduced system is, in general, non-Markovian, and can for a large class 
of processes be described by a generic non-Markovian GME of the form [TJ [HI [22] 
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Here, the reduced density matrix p{n, t) of the system has been resolved with respect to 
the number of transferred charges n. The memory kernel W describes the influence of 
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the environment on the dynamics of the system, while the inhomogeneity 7 accounts for 
initial correlations between system and errvironmentS For a given system, the derivation 
of an equation hke Eq. ([1]) may be a difficult task, but many different examples can 
already be found in the hterature [HI [23l [2l]. Both W and 7 decay with time, usually 
on a comparable timescale. We consider systems where W and 7 with time decay 
faster than any power-law. With this condition it can be shown that the effects of 
the inhomogeneity vanish in the long-time limit. The inhomogeneity 7 is consequently 
irrelevant for all statistical quantities in the long-time limit. For finite times, it may, 
however, play a crucial role [T9] . 

The cumulant generating function S{x,t) corresponding to P{n,t) is defined as 

^s(x,t) ^ J2 P{n, t)e'"^ = Tr{p(n, t)}e*"^ (2) 

n n 

where x is the so-called counting field. The second equality defines P{n, t) as the 
trace of the n-resolved reduced density matrix. The m'th cumulant {{n'^)){t) is directly 
connected to the Taylor coefficients of the CGF in Eq. ([2]) according to the definition 
{{n'^))(t) = d"^S{Xjt) /d{ix)"^\x^o- ^^w derive a general expression for the CGF of 
a system described by a GME of the form given in Eq. ([1]). In Laplace space, defined 
by the transform /(x, z) = J2n Io° ^^fi^^ t)e*"'^~^*, the equation has the algebraic form 

zp{x, z) - p{x, t = 0) = W{x, z)p{x, z) + 7(x, z), (3) 

which can formally be solved for p{x, z) by introducing the resolvent Q{x,z) = 
[z — yV{x, z)]~^ . Returning to the time domain by an inverse Laplace transformation, 
the CGF becomes 

1 na+ioo 

e^(x,t)^ / dzTT{g{x,z)[pix,t = 0)+iix,z)]}e^\ (4) 

where a is a real number, chosen such that all singularities of the integrand are situated 
to the left of the vertical line of integration. This expression constitutes a powerful formal 
result, but as we shall see in the following, it also leads to useful practical schemes. 

As already mentioned above, the CGF scales linearly with time in the long-time 
limit for kernels that decay faster than any power-law and is independent of the initial 
conditions [T8]. For such systems we can define the zero- frequency cumulants of the 
current as ((/™)) = ^((n"^))(t)|t^oo (with e = 1 in the following). With the counting field 
X set to zero, the system tends exponentially to a unique stationary state determined 
by the 1/z pole of the resolvent Q{x = 0, z). The stationary state is given by the 
eigenvector corresponding to the zero-eigenvalue of Wo = W(x = 0,z = 0), i.e., 
limt^ooP(x = 0,t) = |0)), where |0)) is the normalized solution to Wo|0)) = 0. With 
finite values of x, an eigenvalue Xo{Xy z) develops adiabatically from the zero-eigenvalue 
and the long-time behavior is still determined by the pole l/[z — Ao(X)^)] of Q{Xiz) 
close to zero. The particular pole Zq{x) that solves the self-consistency equation 

2^0 - Ao(x,2o) = 0, (5) 

I With environment we refer to all degrees of freedom that have been traced out. 
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and goes to zero with x going to zero, i.e., -2o(0) = 0, consequently determines the long- 
time limit of the CGF. We thus find e*^*^^'*^ T){x)(^^°^^''^ for large t, where 'P(x) is a 
time-independent function depending on the initial conditions. The current cumulants 
then read ((Z"^)) = ^fal^^^ Ix^o- Markovian limit for the kernel W(x, — >■ 0) 

we obtain zq{x) = ^o{Xj^)j consistently with previous results for the Markovian case 
[T6| [T7] . From Eq. ([5]), it is clear that calculations of current cumulants proceed in two 
steps: First, the dominating eigenvalue Xq{x,z) has to be determined. Secondly, the 
self-consistency equation must be solved for zo{x)- 

In the original approach to calculations of CGF for non-Markovian transport 
systems a so-called non-Markovian expansion was developed [18]. In this approach, 
the CGF is expressed as a series using only the Taylor expansion of the dominating 
eigenvalue around z = 0. Mathematically, the non-Markovian expansion is equivalent to 
the solution of the self-consistency equation. The method, however, requires an analytic 
solution for the dominating eigenvalue with its full dependence on the counting field |18j . 
This is not a feasible approach, when the involved matrices are large or, equivalently, 
for systems with many degrees of freedom. However, if we are only interested in a 
finite number of cumulants an alternative route exists. To this end, we have developed 
a scheme for calculating finite orders of current cumulants, where both the eigenvalue 
problem and the self-consistency equation are solved using perturbation theory in the 
counting field. A particular strength of the scheme is that it is recursive, allowing for 
calculations of cumulants of very high orders [19], [21]. Here, we do not present all details 
of the derivation of the recursive scheme, but mainly focus on the final results for the 
first three current cumulants; mean current, noise, and skewness. 

We consider an expansion of the dominating eigenvalue in x and z, Xo{x,z) = 
Srz=o ^Tpfr^'''^''^' ^^^^ suppose that the coefficients c^'^''-' are known. This 

allows us to solve Eq. ([5]) for zo{x) to a given order in x- From the expansion 
-2o(x) = J2'^=i ^^^((-^"))> then extract the zero- frequency cumulants of the current. 
The results for the lowest cumulants are 

an = c^'''\ (6) 

((/2)) = c(^-o) + 2c(^-°)c(i'^), (7) 
= + 3c(2.o)c(i.i) + 3c(^'°) [c(i'°)c(^'2) ^ 2(c(i'i))2 + c(2.i)] . (8) 

Higher-order cumulants are readily calculated in a recursive manner [T9l [21] . The 
cumulants consist of contributions from the purely Markovian quantities c^'^'"-* and 
the non-Markovian terms c*^'^''-* for < /. We observe the general rule that the n'th 
cumulant requires knowledge of non-Markovian terms c^'^'^^ of order < I < n |18j . 
Consequently, the mean current is not sensitive to non-Markovian effects, whereas 
higher-order cumulants are. In general, the distinction made here between Markovian 
and non-Markovian systems is related to the existence of a finite memory in Eq. ([1]) . In 
practice, the memory effects depend on the choice of degrees of freedom that are traced 
out. It may for example be convenient to trace out degrees of freedom of a Markovian 
system and describe the resulting, reduced system by a non-Markovian GME as it was 
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done in Ref. Of course, in such cases our method yields identical results regardless 
of the particular formulation, Markovian or non-Markovian. 

We still need to calculate the expansion coefficients c*-'^''-' entering the expression 
^oix^ = X]r«=o ^rrfr'^'''^''^- '^^^ eigenvalue problem for Xo{x, z) is defined by the 
equation 

Wix,zmx,z))) = [Wo + W'(x,^)]|0(x,^))) = Ao(x,^)|0(x,^))), (9) 

where we have written W(x, z) as the sum of an unperturbed part Wo and the 
perturbation >V'(x, -z) = VV(x, z) — Wq. The right eigenvector |0(x, -z))) reduces 
adiabatically to the stationary state |0)) with x s'^nd z going to zero and satisfies the 
normalization condition ((0|0(x, z))) = 1. Here, the left eigenvector ((0| is defined by the 
relation ((0|>Vo = 0. Using Rayleigh-Schrodinger perturbation theory we can express 
the coefficients in terms of the Taylor coefficients of W'{x^z) denoted as W'-'^''^ i.e., 
^'ix,z) = Em=o Vi^^'''^^ = by definition, and the pseudoinverse of 

the kernel 7^ defined as 7^ = QW^^Q. Even if Wq is singular, the pseudoinverse is in 
fact well-defined since the singular part of the kernel is projected away by the projector 
Q = 1 — |0))((0|1§| We can now calculate the expansion coefficients and we report here 
the resulting expressions for a few of them 

c(i'O) = ((0|W(^'°)|0)), (10) 

c(i'i) = ((o|(w(^'^) - w(^•°)7^w(°•^))|o)), (11) 

c(2'0) = ((0|(W(''°) - 2>V(''°)7^>V(''°))|0)). (12) 

The expansion coefficients can also be calculated in a recursive manner as described in 
Refs. [m [21] . We note that the recursive scheme can be used both for analytic and 
numerical calculations. Evaluation of the pseudoinverse TZ amounts to solving matrix 
equations which is feasible even with very large matrices j25J. Numerically, the scheme 
is stable for very high orders of cumulants (> 20) as we have tested on simple examples. 

3. Dissipative double quantum dot with a QPC charge detector 

We illustrate our method by considering a model of charge transport through a double 
quantum dot embedded in a dissipative environment, see Fig. [1^. A QPC close to the 
double quantum dot is used as a charge detector |26l EZl [28] . The double dot is operated 
in the Coulomb blockade regime close to a charge degeneracy point, where maximally 
a single additional electron is allowed to enter and leave the double quantum dot. The 
Hamiltonian of the full setup is 

H = Hoo + Hqpc + Hdd_qpc + Ht + + Hp + Hb + H^o-b (13) 
where the various terms are defined in the following. 

§ Details of the super-operator notation used here can be found in Ref. [25] . 
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Figure 1. a) Schematics of the double dot and the QPC charge detector. The state 
of the charge qubit modulates the transparency of the quantum point contact via the 
direct Coulomb interaction between them, b) Mean current ((/)) /Tc on a log-scale vs 
level detuning e/T^; line styles correspond to different qubit regimes: coherent case 
{Td = Tr/2, a = 0, solid), coupled only to the QPC (a = 0, Vd/T^ = 0.075, dotted), 
coupled only to the phonon-bath (r^ = ^r/'^, oi = 6.25 x 10^^) for decoupling i) 
(dark grey) and decoupling ii) (light grey), coupled to both QPC and phonon bath for 
decoupling ii) (dot-dashed). Other parameters are T^/Tc — 10, T^/Tc = 10 and, for 
gray and dot-dashed curves only, ksT/Tc = 5. See text for descriptions of decouplings 
i) and ii). 



The Hamiltonian of the double quantum dot is i^DD = f^z + TcSx, introducing the 
pseudo-spin operators = \L){L\ — \R){R\ and s^^ = \L){R\ + \R){L\. The tunnel couphng 
between the two quantum dot levels \L) and \R) is denoted by T^, while e is the energy 
detuning of the two levels. The pseudo-spin system is tunnel-coupled to left (L) and right 
(R) leads via the tunnel-Hamiltonian Ht = Y^j^ a=L Ri^kaCUO){(^\ + h.c), with both 
leads described as non-interacting fermions, i.e., Ha = ^fcacI^Cfc^) = L,R. We 
furthermore include a dissipative environment consisting of a reservoir of non- interacting 
bosons Hb — Ylj ^jO^aj. The heat bath couples to the component of the pseudo-spin 
via the term Hdd-b — VbSz with Vg = J^^- ^(o] -|- where Cj is the electron-phonon 
coupling strength. 

The Hamiltonian of the QPC detector is Hqpc — Ylk a=L R^><:<^^ka^ka + 
X^fefe'-^o dl.^dk'R + h.c, where the first term models the QPC leads and the second 
term describes tunneling between them with a real energy-independent tunneling matrix 
element Tq. The interaction between the QPC and the double quantum dot is given 
by Hqpc-dd = Y.k,k',j=R,L^'^j dli^dk' + h.c. with 6Tj describing the variation of 
the QPC barrier opacity due to a localized electron occupying state j = L,R. If 
5Tj = 0, the current through the QPC at zero temperature is 7^^^ = 2ttT^V lV rc^V / h 
independently of the state of the qubit. Here Pl/r are the density of states of the QPC 
leads and V is the bias across the QPC. If 5Tj ^ and the qubit is in state the 
QPC current at zero temperature is /j^^^ = 27r(To + ST^fVLVne'^V/h. For 5Tr ^ 5Tl 
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the QPC introduces a decoherence mechanism of the charge qubit, because it effectively 
"measures" the right and left states of the qubit. In this paper we consider the weakly 
responding limit where the QPC current is only slightly modified by the charge state of 
the qubit, such that - < I^^^ [29]. 

To describe charge transport through the charge qubit and the QPC we follow the 
scheme outlined in Fig. [1] (see also Ref. [26]). Here, the variable n (m) corresponds 
to the number of electrons counted in the right lead of the qubit (QPC). Instead of 
writing the equation of motion (EOM) for the reduced density matrix resolved with 
respect to the (n, m) components, it is convenient to introduce the counting fields x 
and (p, corresponding to n (qubit) and m (QPC), respectively. We then trace out the 
leads of the qubit and the QPC, leading to an EOM for the reduced density matrix 
= (o"oo, o'LL, d'RR, aiR, ctrlY' of the double dot and the bath of bosons. The elements 
dij are still operators in the Hilbert space of the boson bath. This approach is valid 
to all orders in the tunnel coupling T^. under the assumption of a large bias across the 
system and the QPC charge detector [261 [30]. The Liouville equation is then 

&{t) = £(x, 0)a(t) - i[Hb + Vbs... a{t)l (14) 

where C{x^ 0) describes the reduced dynamics of the charge qubit and the QPC, while 
the other contributions are given by the heat bath Hamiltonian Hb and the interaction 
term Hdd-b = VrSz- The Liouville operator £(x, 0) for the combined qubit-detector 
system is 

/ -Tl + Doh{(j)) Trc'^ ^ \ 

Tl DLh{<f)) iT, -iT, 

-VR + DRh{(i)) -iT, iT, (15) 

iT, -iT, -ie-Vd{(p) 

V -iT, iT^ ie-Vd{(p) j 

with /i(0) = (e*"^ — 1) and the energy-independent rates = 27r ~ ^fcc); 

a = L, R. These rates describe charges entering (leaving) the left (right) quantum dot 
from (to) the left (right) lead. The number of electrons n that have tunneled to the right 
lead is increased by tunnel processes from the right dot with rate Tr, and the counting 
factor e^^ consequently enters the corresponding off-diagonal element of the matrix in 
Eq. (USD. 

The diagonal terms Dj{e^'^ — 1) describe counting of tunneling events in the QPC 
charge detector with the tunneling rate Dj = Ij/e depending on the state of the qubit, 
b); j = 0,L,R. The state |0) corresponds to the double dot without an additional 
electron. The generalized dephasing rate Frf(0) = [TR + {^/Dl — y/DR)'^ — 2\/DLDR{e^'^ — 
l)]/2 represents the decoherence of the off-diagonal terms aLR. and a^L, taking into 
account the dephasing induced by the QPC. For = it yields [F^ + (-\/Dl — v^Ur)^]/2, 
the dephasing rate expected from the coupling of the qubit to the right lead connected 
to the double dot and the nearby QPC charge detector [2S]- If = Dr, the QPC does 
not detect the position of an electron on the double quantum dot. In that case, the 
qubit is not loosing its coherence due to the QPC, but only due to the right lead, which 
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contributes with the decoherence rate The functional dependence on x ^^id 

contains information about correlations in the transport statistics of the combined QPC 
and qubit system. Such correlations will, however, not be explored in further detail in 
this work. Here, we limit ourselves to studies of the transport statistics of the qubit. We 
note that it is also possible to extend the formalism to cases where the thermal energy 
ksT is comparable to the bias V across the QPC p9] . 

We see that the EOM for a{t) in Eq. f|T4|) clearly is a Markovian GME defined on a 
(infinitely) large Hilbert space due to the inclusion of the bosonic heat bath. To reduce 
the dimensionality of the problem we trace out the boson degrees of freedom, and as we 
shall see this leads to a non-Markovian GME. We first consider the electronic occupation 
probabilities pi = Tisio'ii}, i = 0,L,R, where Tr^ is a trace over the bosonic degrees 
of freedom. The dynamics of the occupation probabilities follow from Eq. ( fT4|) 

Po{t) = [-Tl + Doh{^)]poit) + Tne'^pRit), (16) 

PL{t) = TLPoit) + DLh{4>)pL{t) + lTj:iB{(rLR{t) - a RL{t)} , (17) 

PR{t) = [-Tr + DRh{4>)]pR{t) - zr,TrB{<Tiii(t) - <TRL(t)}. (18) 
The effect of the bath enters only via the dynamics of the off-diagonal elements gi^r and 
gri,. The EOM for oj^r can be written oj^r = -X+{(j))d-LR - iiH'^^aLR + olrH^b'') + 
iTc{aLL — o'rr), where A±(0) = ±ie + Td{(f)) and = Hr ± Vr- Formally, this 

equation can be solved as 

<TLR{t)=iT^ /*rfre-^+(<^)(*--)e-^^^"'(*-^)[az.^(r) -a^^(r)]e^^^"'(*-^) 

+ e-"+W*e-*^^^'*aiR(0)e^^^"'*, (19) 

where the term containing the initial condition aL/j(0) eventually enters the 
inhomogeneity [19]. We will only be considering zero- frequency cumulants and can 
thus safely neglect this term. We do not show the similar solution for aRiit), but it is 
important to note that crRiit) is not simply the complex conjugate of aiRit) due to the 
counting fields x ^"^^ 0- Only in the limit 0, the standard relation between the 
off-diagonal elements is reestablished. 

Substituting the solutions for aiRif) and crRLit) into Eqs. f|T7j) and f|T8l) . we can 
obtain a closed system of equations by performing a decoupling of the charge degrees of 
freedom and the boson bath. Two possible decouplings are considered: 

(i) the standard Born factorization, where the system and the bath degrees of freedom 
are factorized as an pi® with ap = e~^^^ / Tr^je"'^^^}. 

(ii) the so-called state dependent Born factorization [25j, where the heat bath is 
assumed to equilibrate corresponding to the given charge state, such that (Jll/r.r — 
pL/R®(^'p~^ ■* with a^iP = e~^^B Tir^c'^^b This is equivalent to the standard 
Born approximation, afler the qubit and heat bath have been decoupled via a 
polaron transformation at = 0. 

Here, (3 = l/Zc^T is the inverse temperature. These decouplings are valid when the 
bath-assisted hopping rates r^^(2;) (proportional to T^) are much smaller than Tl/r. 
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Additionally, approximation (i) is only valid when the strength of the electron-phonon 
coupling is so weak that the state of the qubit does not affect the equilibrium of the 
heat bath a/j. We derive an expression for the memory kernel using assumption (ii). 
The result corresponding to assumption (i) can easily be obtained via the substitution 

The memory kernel 0, -z) for our model, with p = {po, Pl, Pr)'^ , is given in 
Laplace space as 

/ -TL + Doh{^) r^e^^^ \ 

Tl -T^^\z,<P) + DLhi<P) r(5-)(^,0) (20) 

V r(,+)(2,0) _r(j-)(^,0)-r^ + D^/,(0) y 

Most notably, the bath-assisted hopping rates are r^^\z,4>) = T^{g^^''[z±{(j))] + 
with 2;±(0) = 2; — A±(0). The particular dependence on the QPC counting 
field comes about via the dephasing rate Td{(l))- The bath correlation functions g^'^^z) 
in Laplace space follow from the corresponding expressions in time domain which can 
be written as 

g(^)(t) = Tr^{e-^^B^'Vf ^e^^^"'*}. (21) 

Upon the substitution aj^^ a/j, the result for approximation i) is obtained. In that 
case we have r^-'(2;, 0) = r^''(2;, 0). 

We proceed by calculating the bath correlation functions g^^\t) using assumption 
i). It is convenient to change representation in Eq. fl2T|) . writing 

^(±)(t) = TTs{e~'^'B^''e^^-%e-'^-'e'^B''} 

= TrB{f//^)(t)a^f/(^)(t)}, (22) 

where we have introduced the operators U{±){t) = e'^^^^e^^s and moreover used 
the fact that ap does not evolve with the stationary Hamiltonian Hb- It is easy 
to demonstrate that the EOM for the operators U{±){t) in this representation is 
5if/(±)(t) = ^iVg(t)?7(±)(t) with VBit) = e-'^B^Vse'^Bt^ ^^n thus identify these 
operators with the evolution operators in the interaction picture. The solution of the 
time dependent differential equations is f/(±)(t) = T[e^*-^ '''^^s*^'^^] with T[-] being the 
time-ordering operator. Equation ( |22l) can then be written 

g^^\t) = Trij{T[e^'^i'o'^-^sM] f\^^Ti S^, drV^ir)^-^ 

= TrB{Tx[e^'^o'^*'"^^^(*'^]}, (23) 

where we have introduced the time-antiordering operator T[-\ and the Keldysh contour 
ordering operator Tft:[-]. The second equality expresses the bath correlation function 
in terms of a Keldysh propagator with a branch dependent interaction potential 
T^VB{t') with r3 = +(— ) for the forward (backward) Keldysh branch. This expression 
can be evaluated using perturbation theory in the electron-phonon couplings Cj and 
we obtain g'^if) = 1 — dojJ{uj){[l — cos(co')f:)] coth(/3tt;/2)}/co'^ + 0{\cj\'^), where 
J{oj) = J2j |ciP<5('-^ — ^^j) is the spectral density of the heat bath. Below, we consider 
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the case of Ohmic dissipation J{uj) = 2auoe We note that g''~^\t) = g^ \t) = g(t). 

The Laplace transform of g{t), to leading order in 1/ (3uJci reads 



g(z) = -^{l-2a 



(24) 



where is the digamma function and T>{x) = ci(x) cos(x) + sin(x)[si(x) — 7r/2] 

with the sine and cosine integrals defined as si(x) = dt sm{t) / 1 and ci(x) = 
— dt cos(t)/t, respectively. 

For decoupling ii) we can calculate the bath correlation functions exactly to all 
orders in a using standard many-body techniques ^21j. The results for the bath 
correlations functions are g^^\t) = e~^^^*\ where 

W{t) = dcu^ I [1 - cos(cut)] coth (^^^ + % sin(cut) | . (25) 

In the weak coupling limit, we can expand this expression to first order in a. The 
bath correlation functions cf'^^{z) for an Ohmic spectral density, to leading order 
in l//3u;c, are g'^'^'^^z) = g{z) ± i2aF{z/uJc)/ z where g{z) is given in Eq. (j24l) and 
F{x) = ci(x) sin(x) + cos(x)[7r/2 — si(x)]. We see that the difference between the bath 
correlation functions obtained within the two decoupling schemes is proportional to the 
imaginary of part Eq. (l25l) . This difference is the main reason that decoupling i) does 
not lead to any asymmetry of the cumulants between the emission and absorption sides, 
as we shall see. We note that the bath correlation function is g{z) = 1/z for a = 0. 
This corresponds to the coherent regime of the qubit and the 2r-dependence of T^^\z, (j)) 
describes the "effective" memory due to the off-diagonal elements that have been traced 
out. Furthermore, for ^ ^ the rates become rg'^(O,0) = 4:T'^Vd{(t))/[Td{(t)f + 45^], 
and for = we obtain the standard expressions for incoherent tunneling rates in a 
double dot 

We now consider the current cumulants of the charge qubit for = 0, where the 
effects of the QPC is captured by the dephasing rate F^. Evaluating Eqs. ([6]) and (IH]) 
using the kernel in Eq. ( l20i) we find the analytic expressions for the cumulants. The 
expressions for the current and noise with a = coincide with known results for the 
coherent case [301 [3l], however, with the total dephasing rate being F^ = + {y/D^ — 
V-D/j)^]/2 rather than just F/j/2. Considering only the Markovian contributions c^^'^\ 
we obtain the well-known analytic expressions for sequential tunneling [[32]. We next 
compare different regimes with various strengths of the QPC dephasing rate F^ and 
different electron-phonon couplings a. In Fig. [lb, we show the mean current as function 
of the level detuning e. The solid black line corresponds to the coherent regime, where 
the peak is symmetric around e = and the width of the peak is proportional to F/j/2. 
The dotted hne shows the effect of the dephasing due to the QPC charge detector, 
without any contributions from the boson bath (a = 0). For ^ Di, the peak width 
is given by the total dephasing rate, i.e., the sum of the intrinsic contribution F/j/2 and 
the contribution from the QPC [y/Di — a/Dr)^/2. The width is thus larger than in 
the coherent case. The peak, however, remains symmetric. The gray curves show the 



The influence of charge detection on counting statistics 12 




Figure 2. a) Fano factor ((/^))/((/)) as a function of the level detuning e/Tc- b) 
Normalized skewness ((/■^))/((/)) vs e/Tc- Parameters corresponding to the different 
line styles are given in the caption of Fig. [1] 

mean current, when the double dot is coupled to the heat bath at finite temperature, 
but with the QPC charge detector not detecting the position of an electron on the 
double dot {Dr = Dl). The two gray lines correspond to the two different decouplings, 
decoupling i) with light gray and decoupling ii) with dark gray. We note that the curve 
corresponding to decoupling ii) is asymmetric around e = 0. Finally, the dot-dashed 
curve corresponds to the mean current with contributions to the dephasing from both 
the QPC and the heat bath. 

While the mean current studied so far only reveals little information about the 
dephasing mechanisms, we expect more information to be contained in the higher-order 
cumulants. In Fig. [2] we show the Fano factor ((/^)) / ((/)) and the normalized skewness 
((/^)) / ((/)) as functions of the level detuning e. Line styles and parameters are the same 
as in Fig. [T]d. Generally, we observe that the two cumulants to a higher degree than the 
mean current discriminate between various regimes. In particular, the coherent regime 
has strongly super- Poissonian behavior with both of the normalized cumulants reaching 
values larger than the Poissonian limit of 1. In contrast, for the sequential tunneling 
regime, given by the Markovian contributions to the cumulants, only sub-Poissonian 
behavior is observed [32] . 

If a dephasing mechanism is introduced, either due to the QPC charge detector or 
the heat bath, the super- Poissonian behavior is gradually reduced towards the sequential 
tunneling regime with increasing dephasing rate. The dephasing rate P^ due to the QPC 
can be modified via the coupling between the charge qubit and the QPC. On the other 
hand, the dephasing induced by the heat bath changes with the bath temperature, 
which can also induce the transition between coherent and incoherent transport. Such 
a transition was recently observed in shot noise measurements of transport through a 
double dot at low temperatures [32]. Comparing the scales of the vertical axes of the 
second and third cumulants in Fig. [21 we conjecture that such a transition may be better 
visible in the third cumulant and that even higher order cumulants in general could be 
more sensitive to such a transition. 
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We now study, how the symmetry of the cumulants changes in the various regimes. 
Without couphng to the QPC or the heat bath the cumulant are symmetric around 
e = 0. In case the qubit is coupled only to the QPC, this symmetry remains intact. 
However, with non-zero coupling to the heat bath, asymmetry around e = is observed 
(see dark gray and dot-dashed lines). The asymmetry occurs due to the asymmetry in 
emission and absorption of bosons at low temperatures. Phonon absorption and emission 
dominate for £ < and e > 0, respectively. We note that the curves obtained using 
decoupling i) do not capture this essential physics, not even to first order in a. Indeed, 
in all the regimes for which the cumulants are symmetric, the rates have the property 
r*^+^(2;) = T^~\z). This particular property is not fulfilled in decoupling ii) due to the 
imaginary part of Eq. fl25|) . and we believe that decoupling ii) should be correct to any 
order in a. What is then the essential difference between the dephasing induced by the 
QPC charge detector and the bath-induced dephasing in the model we are considering 
here? The QPC charge detector induces dephasing with little influence on the dynamics 
of the system and does not destroy the symmetry in the coherent regime. In contrast, 
the bath-induced dephasing is generated with emission and absorption of bosons, and 
the charge qubit is influenced by the intrinsic asymmetry of this process. We conclude 
the analysis of this model for now and postpone a further analysis to future research 
aimed at an improved understanding of these dephasing mechanisms. 

4. Outlook and open questions 

We have shown how it is now possible to calculate cumulants, also of high orders, for a 
wide class of counting problems. In this work we have applied the method to investigate 
how the dephasing mechanism of a QPC charge detector differs from that of an external 
heat bath. The method is, however, applicable to a large class of non-Markovian 
counting problems, also from outside the field of physics. The general approach outlined 
here is also suitable for calculations of the finite frequency noise spectrum [19]. It 
is relevant to ask if it is possible to extend this approach to frequency dependent 
cumulants of higher orders, similar to the results already obtained for Markovian systems 
[33] ? Another interesting issue concerns systems with very long memory time (for 
example. Levy- flight processes). Such cases would require us to reconsider some of the 
assumptions underlying the theory presented here. 

We would like to thank R. Aguado, T. Brandes, A. -P. Jauho, S. Kohler, K. Netocny, 
and M. Sassetti for fruitful discussions and suggestions. The work was supported by 
INFM-CNR via "Seed" project, the Villum Kann Rasmussen Foundation, the grant 
202/07/ J051 of the Czech Science Foundation, and research plan MSM 0021620834 
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